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ABSTRACT 

This paper addresses the breakdown of the classical (CBES) field aligned 
transport relations for electrons in an inhomogeneous, fully ionized plasma as a 
mathematical issue of radius of convergence, defines for the first time the finite 
Knudsen number conditions when CBES results are accurate, and presents a 
global-local (GL) way to describe the results of Coulomb physics moderated 
conduction that is more nearly appropriate for astrophysical plasmas. This paper 
shows the relationship to and points of departure of the present work from the 
CBES approach. It also presents an analytical, more tractable, and economical way 
of evaluating the global aspects of the collisional physics than has been done 
numerically in the previous paper . The CBES heat law in current use is then shown 
to be an especially restrictive special case of the new, more general GL result. 

A preliminary evaluation of the dimensionless heat function discussed in the 
previous paper (Olbert, (1983)), but using analytic formulas, as presented below, 
shows that the dimensionless heat function profiles versus density of the type 
r' - — — aw cm or! Vii dh 9HPpH Sdl 3 T W1 fid 300631* possible . 

Introduction 

It has long been recognized that the Chapman ( 19 l 6 )-Chapman and Cowling 
( i 939 ) -Brag inski i ( 1965 )-Enskog ( 1 917 )-Spit zer-Hftrm ( 1953)-Spit zer (1962) (CBES) 
transport results become unreliable when the Knudsen number, K, (a kinetic free 
path over macroscopic scale length) is not small. A serious outstanding issue in 
this regard is how small must K be in order that these transport formulas are 
accurate. This difficulty plagues the theoretical investigation of conduction 
supported, fully ionized stellar winds. The usual interpretation of this 
breakdown is that relevant physics has been omitted. Wave particle effects are 
often suggested to be important in these locales. However, this conclusion is not 
inescapable. Subtle prejudices have, until now, led to the insistence that the 
physics of these regions must remain local, even if wave particle effects must be 
invoked to do this. These assumptions are largely shaped by our everyday 
experience and are fortified by attendant mathematical convenience. The physics 
of Coulomb moderated transport is not amenable to this _ad hoc approach as we have 
previously shown (Scudder and Olbert 1979a, b, SO I, II, Olbert ( 1 983 ) ) • 

Since electrons carry almost all the heat flux in a plasma the kinetic 
discussion that follows is oriented toward electron transport signatures. Because 
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the bulk flows of astrophysics are usually smaller than characteristic electron 
thermal speeds certain dimensionless ratios are tailored for this electron 
situation. Unless specifically stated otherwise we will refer to electron kinetic 
properties for the remainder with the major exception of the plasma conservation 
equations where certain results pertain to the fluid as a whole. Because 
astrophysical plasmas are strongly "magnetized” we confine our attention to field 
aligned transport signatures. 

1 • The CBES Postulates 


The CBES transport description for electrons proceeds under the interlocking 
postulates 


1) that the microscopic electron velocity distribution in the proper frame, 
f e (w). that solves the kinetic equation is separable into a sum of two terms: 

CBES 

f e = f e +6f e (K) (1) 

such that all finite free path effects are solely contained in the Knudsen number, 
K, functional dependence of 6 f ; 

2) that the velocity dependence of f^ ^ is only through proper frame speed w = 

I w| ; 

3) that f° is a local Maxwellian distribution: 



(where k is Boltzman f s constant, n, T, and m are the local density, 
temperature and mass, respectively); 

4) that 6f is a perturbation of f ^ such that 
e e 


6 f (w) 

* 


0(K) , 


for a sufficiently broad range of speeds to control the transport 
signature sought; and 


(3) 


5) that the perturbation ordering parameter is the Knudsen number, 

K, which is assumed infinitesimal and defined as the ratio of a kinetic 

scale, X , to a macroscopic scale, L M . The macroscopic length L,., defined 
by MM 


L 


-1 

M 


d 

| — in M| , 
ds 


(4) 


is the shortest scale length along the arc length s of the magnetic field among 
the macroscopic variables, M: density, n, temperature, T, magnetic field 
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strength, B, or bulk speed, U. In the astrophysical 

the density profile; we assume this in the remainder. 

strongly speed dependent in a plasma. It is usually chosen as that spatial length 

within which substantial isotropization and approach to a Maxwellian takes place. 

Those two physical time/length scales are commensurate (cf. Spitzer, 1962, p. 1 33 ) 

for particles with speeds w_ equal that 

1 *6 


context L m is usually set by 
The kinetic scale, X, is 


of the root mean 


square speed 



Spitzer (1962) defines the related 


e 


characteristic time, t , (eqtn 5-26) for this process for electrons to yield a 
Spitzer kinetic scale Free path of the root mean square (rms) electron, 


X rms 5 w T,e '“T.e 1 

- 

8(0.71*4) tne 4 Hn A 


where In A is the Coulomb logarithm. Braginskii (1965) defines "the 
electron... collision time" with a formula (cf. 2.5e) that differs from 
Spitzer' s by 3%- According to this practice "the Knudsen number" defined by 


K = 


X rms (w T,e ) 


( 6 ) 


L 

n 

appears as the dimensionless parameter for the plasma style Chapman-Enskog 
expansion . 


The CBES approach attempts a perturbation expansion for 6 f of the form 


6 f 


CBES 


■ <*<V L Te- 


L n /L U’ 


L n /L B’ 


U, w) K + 0 


(K 2 ) 


(7) 


aria tne approximation is usually truncated with the first term in Cf). Note tnat 
the structure of ^ depends only on the velocity space w, on local macroscopic 
.variables and ratios L /L M of the scale length of the density to those of other 
macroscopic variables evaluated at the point x. Consequently 3* i- s independent of 


K. 


The CBES truncation (7) of 6 f translates immediately into a similar 
truncation of a Maclaurin series For Q„ since 

Q„ = /// (f (0) + 6f ) 1/2 m w 5 cose sine d0d4>dw. (8a) 

" e e e 

By postulate (2) f ^ is even in velocity space so that (8a) reduces to 

Q„ = 0 + Iff (6f g ) 1/2 m g w 5 cose sine ded 4 )dw. (8b) 

Using equation (7) we obtain 
CBES 

Q„ = K Iff 3 1/2 m w 5 cose sine ded$dw. (8c) 

The remaining CBES postulates allow a reduction of (8c) to functional 
dependence on lower order moments; this permits a heat law for the fully 
ionized plasma, parallel to Fourier's law, that permits closure, viz: 
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(9a) 


CBES 

Q„ “ 0 + 


% T e 5/2 


dT 

e 

ds 


which is easily rewritten as the infinitesimal heat flux 


CBES V T 7/2 L 
Q„ « — ] K, 

x rms L Te 

where for a proton-electron plasma (Spitzer (1962)) 


(9b) 


,7/2 


X = (0.419) (0.225) 20 (-) 3/2 


( 10 ) 


1/2 4 

ro e e in a 


Notice that the flow of heat in the CBES approach (8c) arises out of the non-zero 
perturbative correction to adiabatic (Q n = 0) behavior. 


If the assumptions made to derive (9) are violated this marks the failure of a 
mathematical artifice: a simple, straightforward Maclaurin series perturbation 
scheme. It is, for example, well known that the Chapman-Enskog scheme does not 
provide access (even in principle) to the broadest class of solutions to the 
Boltzmann equation (Grad (1958)). Coulomb moderated transport effects are 
unavoidably present, but inaccessible within the CBES mathematical machinery. 
Extensions of the CBES approximation after the form of Burnett and others are 
analogous to attempts to enhance the radius of convergence of the Maclaurin series 
for the velocity distribution function, f(w) , and, therefore Q„ beyond those in 
equation (9) by retaining additional terms in the power series. These procedures 
, are cumbersome. This is made clear by looking at a mathematical model of the 
, perturbation issues at hand. 


2. Model Perturbation Issues and their Resolution 


Consider a system whose physical behavior is exactly described by the 
following function 


g(K) = tanh(K) . K>0 (11) 

[In the heat transport problem the exact g(K) a Q„ (K) is an unknown solution of an 
integro-partial differential equation]. If a description of the physical system 
were desired for the K<<1 regime a one term Maclaurin series expansion for tanh(K) 
would yield nearly indistinguishable results for g(K), viz 

g(K) + g (1) (K) = K K<<1 (12) 

/ 1 \ CBES 

(In our model problem g iK) is proportional to Q„ .) If equation (11) (the 
solution of the integro-differential equation) is the complete and accurate 

description of the physical system even for large K, the fact that g^^(K) is 
inaccurate as K gets large implies nothing about the completeness of the physics 
embodied in equation (11). 
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To parallel the multi-term expansion of f (w) , and, therefore Q„ , in the 

higher moment methods, consider additional truncated Taylor expansions g^ n ^(K) for 
g(K) 


g 

g 

g 

g 


( 1 ) 

( 2 ) 

(3) 

(4) 


(K) 

(K) 

(K) 

(K) 


= K - K V3. 5 

= K - K f/3 + 2K c/15, 7 

= K - K J /3 + 2K V15 - 17K '/315. 


(13) 


Any finite truncation g of g chosen in (13) is a polynomial in K which must 
diverge for large K in contradiction to the boundedness of g(K) as K + ® . In any 
case the Maclaurin series for g(K) of course diverges for K >ir/2. This. result is 
known because the full Maclaurin series of g(K) is in hand. “Clearly g^ a is a 
desperate way to approximate g - even for K < ir/2. In constructing perturbation 
solutions that approximate Q„ the full expansion analogous to ( 1 3 ) is not 
available. The finite radius of applicability of the CBES expression for Q lt needs 
to be defined as we do below for the first time. It is quite common in 
non-linear , physical problems that boundedness of certain quantities such as total 
energy or the magnitude of the heat flow regardless of the size of K available 
impacts the forms of perturbation that are allowed and drive the physicist to seek 


more tractable approximations such as g^ provides for g regardless of K 
described below. As such, this feature of this model is realistic. 


CBES 

The inapplicability of Q fl in the finite Knudsen number regime is 
essentially a failure of the above type: namely of a straightforward, local, one 

term, Maclaurin series expansion in K. (Recall that elementary discussions of 
heat flow e.g., Reif (1965) give intuitive derivations based on the energy flux 
imbalance across a mean free path element of the medium centered on the observer.) 


Hi f f Ifj M r .4- t i. « . .. _• / n 1-N . , , 

w v* wmi VUVJ.VU uApuuuxuiio \ w .1 • , u a iv ci V I y l D J J die 

possible that attempt to bridge asymptotic regimes in a way that truncated Taylor 
series cannot. An example of this for the model problem is to use a Pade 
approximant (which is by no means unique!) for g(K) given by 


g P (K) 


K + cK 2 + dK 3 
1 + aK + bK 2 + dK 3 ’ 


K > 0, 


(14) 


with a, b, c, and d, chosen to optimize the approximation for intermediate K such 
as by matching the Maclaurin series of (14) to that of g as far as the free 
constants allow. Notice that for small K, i.e., K « 1/a 

g P (K) * K K << 1 

but that for K » (1/d) 1/2 

lim g P (K) = lim g(K) = 1 

K -*• oo K -*• oo 

p 

Also note that g (K) is not a polynomial in K. In a practical sense (14) 
represents a physically viable, if approximate, re-summed form of the infinite 
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Maclaurin series for g(K) which is well behaved in either the small or large K 
regime as well as for K > ir/2, the radius of convergence of the Maclaurin series 
for g. Obviously, the exact solution (11) is preferable to (14), if available. 

In most realistic physical problems the exact solution is not available in closed 
form. This implies that physical arguments as for example the boundedness of the 
heat flux even for large K can guide the ansatz. 

This paper reports on an initial attempt to provide the analogue of equation 
(14) for Coulomb physics controlled heat flow Q„(K) for finite K, whose leading 

CBES 

order Maclaurin series reduces for K <<< 1 to Q„ 

3. Established Kinetic Facts vs. CBES Postulates 

At the Solar Wind 4 Meeting in 1978 at Burghausen, we reported on the 
properties of velocity distributions shaped by the Coulomb cross section without 
making gradient expansions (cf. SO I, II; a portion of the previous paper in this 
monograph reports on extensions of that work) . We have demonstrated that 
non-local aspects of the steady state distributions of material along the tubes of 
force play a vital role in shaping the speed dependence in the trans , (E > kT Q ) , 
and extra-thermal, (E > 7kT,) , electrons; the latter group, by observations, 

controls the magnitude of the net skew of and, therefore, the magnitude 

of the heat flux (eqtn 8a, b). 

In the presence of general finite gradients, Maclaurin series perturbation 

expansions for f °^ S (w) about a local Maxwellian are doomed because the strong 
speed dependenceof the Coulomb cross section allows a global communication for 
electrons E > 7kT along the magnetic tube of force (cf. SO I). The information 
that f (corona) >> f (observer) manifests itself in the over populated, skewed, 
tails of the local diltribution which observationally exceed a local Maxwellian 
population by several orders of magnitude above 100 eV at 1 AU. cf. SO II Figure 
1. These are Knudsen number effects of the system's distribution of material in 

the even part of f The speed dependence of f has been observationally 

and theoretically shown to be non-Maxwellian as a result of Knudsen number 

effects. This violates CBES postulate 1 and 3- In this important trans and 

extrathermal regime 6 f/f., ... is not even the same order of K local, which at 

° Maxwellian 

1 AU for example is ^ 1/2, much less infinitesimal as assumed by CBES. 

4) Conservation Laws vs. CBES Implications 

Provided Q„ is not zero, conservation laws for the plasma as a whole 
neglecting viscosity and ohmic entropy production (cf. eg., Rossi and Olbert , 

1970) in steady, field aligned flow imply that 


(3PV)/(2Q„) — in (P/ p 5/3 ) :--tnW. (15) 

ds ds 

In the above P is the total gas pressure of electron and ions, P = P + IP. , p is 
the total mass density, V the magnitude of the field aligned flow velocity, Q„ the 
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I 


I 


I 


f 


I 


parallel heat flux, and W is the dimensionless heat function defined by 

w = Q„/(pVU 2 ) = Q„/(aBU 2 ), (16) 

CO CO 

where U is the asymptotic flow speed and a is the flux tube stream constant 
defined"by the ratio of mass to magnetic flux: 

pV= a §. (17) 

As discussed by Olbert (1981) using current notation (and others) in steady 
state there exists an energy streamline constant defined by 


' o " 

U = - V + + 4» + — , 

2 2 p pV 

where i|> is the gravitational potential. Dividing this expression by the 
square of the asymptotic radial wind speed U and using (17) yields an 
equivalent streamline constant where the dimensionless heat function is 
involved : 


2 1 2 5 P * 

(u/ir) = - (v/u r + ^ ^ + w. 

co oo y s' 

2 2 p U U 

OO 00 

At 1 AU the leading term is the largest and dominant term. At the base of the 
corona the last term is the dominant term. The viability of the coronal 
expansion via conduction reduces to showing that 

1 V(1 AU) 2 

- ( ) ^ W (1.03 R ) 

kP o 

2 U 

CD 

In this sense the size of the dimensionless heat function W is crucial to the 
energetic feasibility of a conduction supported solar wind expansion (cf. 
discussion below). 

Equation (15) may be rewritten without further approximations as 

( PV/Q„ ) = (2 L d )/(3 L w ) (18) 

where the spatial scale for departures from adiabatic expansion, the "diabatic" 
scale, is implied by 


L 1 = | Hn (p 5/3 /P)| (19) 

ds 

and the spatial scale for heat divergence, is implied by 

i d 

L„ =1 In W| (20) 

w . 

ds 

In the CBES perturbation scheme, gradients of all macroscopic variables, M, 
are assumed to be weak in the sense that 
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( 21 ) 


= 0(K) 

rms M 


is sufficiently small to ensure the validity of all the postulates. For a given 
experimental regime of K the practical question arises if it is small enough to 
assure quantitative precision of the CBES formulas. We proceed now to delineate 
what implied conditions on other macroscopic variables are sufficient to guarantee 
the validity of the CBES results. 


Since D and W are combinations of macroscopic variables, CBES ordering (21) 


implies that the associated inverse scales, and are also 0(K) in the 

sense of (21). Within the CBES framework the general conservation of energy, 
(18), can only be consistently satisfied if the enthalpy and parallel heat flux 
are of the same order: 


pv cbes / 




CBES 


= 0(1) 


(22a) 


This constitutes a restriction on the physical problems for which CBES is 
internally consistent since (22a) together with (9b) imply that the enthalpy flux 
must be proportional to K for small K: 


PV 


CBES 


K 6 


(22b) 


Except for the consequences of the CBES postulates (21) there is no general 
ab initio argument that all physical systems will satisfy conservation of energy 
in such a way that the enthalpy flux is comparable to or less than 

CBES 

the CBES heat flux in the Knudsen regime (K «< 1) where Q„ calculation is 

CBES . 

rigorous. By explicitly determining & and, therefore, the limits on V in 

CBES 

the CBES regime, we argue that V so i ar wind ^ ^ throughout the entire solar 
wind expansion. 

It is worth noting that CBES transport is usually discussed (cf. for example 
Cohen et al . (1950)) for plasmas at rest so that the speed, V, vanishes 
identically. Correspondingly (22) is satisfied in the CBES approach as a result 
of a postulate (V i 0) which is generally not true in astrophysical contexts. For 
the remainder we proceed on the assumption that V does not vanish throughout the 
system . 


The CBES heat flux may be rewritten as 


Q 


It 


P 

e 


w K, 

o 


(23) 


where 


w 

o 


3-3 


W 


w 


T.e 


(24) 


a flowing medium to physical 
systems such that the inverse diabatic scale, . is second order small in 


It then follows that the CBES regime refers for 

liabat 

K, (Rossi and 


formal expansion parameter, 
(18, 21-24) imply 


Olbert , 


1970, p. 436) , 


the 

since equations 
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CBES 


where 

In the above the 
given by 


X rms /L D S G K • 


( 25 a) 


G = 1,1 (2P e /P) (T p,n /T e,T ^ (L n /L W J * (25b) 

e 

rras electron transit time of its temperature scale, T is 

e,T e 


%, * * L T/“T,e. 

T e 

e 


and the fluid transit time of the electron density scale height, t 
given by 


is 


t = L /V. 
P.n n 


For naturally occurring systems L is usually larger than L ; however w^ 

e n T.e 

is usually larger than V and the ratio of the two times is more nearly 
commensurate than either the ratio of scales or speeds separately. On this basis 
G is not an enormous nunber being clearly £ 1 for CBES ordering. Using the solar 
corona as a prototype for evaluating G we obtain G < 1/6 beyond the sonic point, 
where we have used the relationship that dfcnT /ddn n £ 1/6 (Sittler and Sc udder 
(1980)) and that dtn W/dtn n £ 1/50 (Olbert (?983) ) . e These considerations 
reinforce the important conclusion that indeed 


X rms /L D 


CBES 
'< 


K 


( 26 ) 


which is a more restrictive, but, nevertheless, consistent ordering of scales 
with the initial CBES ordering (21). 

Apart from trivial constants, the inverse scale length, L f o _1 , associated 

with the change of the electron velocity distribution function at zero proper 
frame velocity is the same as that of the diabatic expansion (26). Therefore, 
(26) may be restated in microscopic language to establish contact with SO I, 
II: the CBES description is internally consistent with conservation of energy 

for systems (such as thermodynamic equilibrium (!)) with L f o that is much 

larger than any of the driving gradient scales of the macroscopic variables, 
l m * Rewriting in symbols equation (26) the inverse scales of f° must be one 
order weaker than those of the macroscopic variables. 

L “ 1 v; (L M _1 ) K. 

This is precisely what is not realized in astrophysical plasmas and leads, via 
the Coulomb cross section "window”, to global lowest order signatures in the 
distribution functions and correspondingly in the transport SO I, II. 
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5 . Enforcing CBES Postulate (3) in an Inhomogeneous Plasma 


In order that the lowest order f be Maxwellian in an inhomogeneous, fully 
ionized plasma (CBES postulate 3), it e is no longer sufficient to only require that 
the rms free path, X , be small compared to the scale height. When K is finite 
in a fully ionized plafma then the rms free path of the thermal and suprathermal 
electrons are no longer infinitesimal, being related to the rms free path by the 
formula 

X (u) = u 4 X , (26a) 

rms 

where u = w/w_ . By contrast collisions between neutral m olecule s have a free 
path with essertfially no speed dependence for speeds above /3kT/mi Meyer (1899)* 

If, as is usual (cf. Parker (1964)), one assumes K << 1 is sufficient for 
CBES heat transport results to ensue then (26a) implies that there always exists a 
finite Uq so that the following ordering of scales is possible: 

X «< L < X( Uri ). (26b) 

rms n (j 

Unlike the situation for collisions between neutrals, equation (26a) implies that 
no matter what the size of the Knudsen number (so long as it is finite) , there 
always exists a speed domain ( u > u ) of the velocity space for which the 
(shortest) macroscopic scale length L is shorter than the free path of members of 
the velocity distribution which have speeds in this domain . This domain need not 
be the sole province of the cosmic rays! 


Having established (26b) it is important to assess the potential impact it 
may have on the kinetic determination of transport in the system. If K be so 
small that u Q is a relativistic cosmic ray and the mean energy of the gas 
is ultra non-relativistic , the practical impact of (26b) on the CBES regime is not 

severe. However, u needs to be well above the dominant proper frame speeds where 
/ Q 

'the perturbations to f ^ kinetically determine the heat flux (equation 8 ), the 
highest moment in the Snsager relations. The precise size of u^ so as to assure 


that non-local effects are collisionally destroyed and do not dominate the 
transport becomes the prerequisite for allowing postulate 3 of CBES to be 
reasonable: find u^ such that 


X 

rms 


< x <y « L „- 


(26c) 


Based on the theoretical properties of the Coulomb cross section SO II argued 
that significant non-local suprathermal particles will be present above the 
dimensionless speed u = 2.16. The observed solar wind electron velocity 
distribution function was shown to develop strongly overpopulated skewed 
suprathermal tails (the halo) above this speed. As an empirical fact the 
in situ heat flux is determined by the skew of the electron distribution near but 
above this transition. We suggest that 


u Q > 3.5 (266) 

is appropriate for the inequality in (26c) . Using this result in £2^c) and 
noting (26a), we can restate the sufficient local condition that f be locally 
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. inhom ° genous * fully ionized plasma without postulating special 
spatial distributions of the plasma in the adjoining volume, viz 8 P 


150 K 


CBES 


« 1 . 


(27) 


otice that this result is considerably more restrictive than the usually stated 
(and implied sufficient) condition in the astrophysical context for CBES validity: 
of < CBES E<?U ^ 10n (27) deflnes the regime for finite K that is consistent with use 


. Froi gB!§ e . infcerna l energy equation, (25), in the consistent 
regime K is algebraically equal to 


Knudsen number 


CBES d , 

K = - (L w /L d ) (P/P g ) (V/w q ); 

(28) 

however, as' the prerequisite for fulfilling CBES postulate (3) 
must be satisfied. This yields the compatibility condition 

equation (27) 

. CBES 
Mg «< 1 * 1 ’ 

(29) 

where the adiabatic Mach number is defined as 

= V/(5/3 P/p) 1/2 

and 

(29b) 

^7 = 0.73 (2P /P) 1/2 (L/L u ) 
1 e D W 

(29c) 

which follows from (28) using (27): 


CBES 

K «< 6.6 x 10“ 3 

(30) 


Equations (29) and (30) are not physical constraints but are coroll aries of the 
ma thematica l postulates made for the CBES approach whirh r»nHo- (-i- ?t ron) 
transport problem to be "local” • ’ “ 


. . ? n tu SS 3 fully ^ on f ze d plasma is prepared in a contrived way (as described 
below) the contradiction of either (29, 30) implies that CBES is no longer 
ppropriate for electron signatures. Depending on the observational technique(s) 
one or the other of these constraints may be more accessible for a distant 
observer of astrophysical phenomena. 


For specially prepared equilibria (as discussed by CBES) such as an 
isothermal flux tube that has slight temperature and density variation but still 
superad iabatic" in the sense of (26) and, therefore possessing very low 


contrast in f (o) between its extremities (a very difficult idealization to achieve 

an astrophysical setting) condition 30 can be relaxed by several orders of 
magnitude to 


isothermal (s) « 1 ’ 

p rovided this cond ition is met everywhere along the magnetic tube of force 
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connected to the point of interest . This is the system to which CBES refers. The 
relaxation of (30) results from a very special preparation of the inhomogeneous 
system. For this circumstance individual trans and extrathermal particles still 
have a global range over arc lengths large compared to * rms ; however, the 
probability function f° refers to the pattern of phase space occupation and it is 
not modified in lowest order in the presence of finite K<<1 provided the system is 
prepared in this special way. Alternatively, in the presence of an absolutely 
uniform T and a macroscopic density gradient the individual trans and 
extrathermal electrons give up sufficient .energy to the polarization electric 
field that the lowest order pattern for i can still be self-similar along B 
even in the presence of strong density gradients. 

We submit that these special systems are so contrived that they are unlikely 
prototypes of plasmas in the astrophysical or other contexts and that conditions 
(29, 30) are more generally applicable. Gravity destroys the weak density 
gradient regimes; T can not be artificially assumed constant throughout a given 
astrophysical flux tube and the lowest order f(w) becomes global in character . 

For the solar wind it is clear by direct observation that at 1 AU that 
Mg >> 1, that K ^ 1/2, that both the inequalities (29 and 30) are strongly 
violated, and that the observed f (w) is not deformed as predicted by CBES theory 
(Montgomery, 1972). What is forcefully driven home by (30) is that the CBES 
description of thermal conduction, for example, is invalid at the very base of the 
coronal expansion. Given observed coronal velocities in the range of 10's km/s 

and an implied base sonic Mach number at 10^° K of approximately .08 and the 

determination in the previous paper that L << L in approximately the ratio of 

1:50 implies (if P = P.) that n W 

e 1 

7 Y[ Q = 5 X 10- 3 

and that 

Mo .r .08 > 

and that the CBES determination of supply of conduction flux to the expansion is 
in error from the base of the corona . This conclusion should be contrasted with 
Parker's (1964, eqtn 76) original argument that CBES should be appropriate for the 
initial portion of expansion possibly as far as 1 AU on the basis of an argument 
like 



"CBES Flux" 



T 5/2 
e 


dT 3 

— - <« - NkT 

e 

dr 2 


3k T 

( -) 

m 


1/2 


"Saturation Flux" 


which is equivalent (using (6, 9b, 10, 24) to 


^Parker 


« .45 (— ) x 2.5 
L N 


(32) 


which we now recognize as an overly generous upper bound for K when contrasted 
with (27). 
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Use of the CBES heat transport description within the conservation equations 
for a flowing medium with M<, > << 1 (eqtn 29) is clearly inconsistent with the 

entire set of postulates necessary to derive equation (9). The failure of such a 
kinetically inconsistent model to predict the observed solar wind properties 
implies nothing about the importance of conduction in driving winds. This point 
is especially crucial since the properties of the conduction flux near the 
supersonic transition are extremely critical for the determination of the 

essentially asymptotic wind properties that are accessible to spacecraft 
observations. 

6. The Global-Local Formulation of the Transport Problems 

In the present restricted monograph we only outline the philosophy and 

present the results of the procedure used to obtain the heat flux form quoted 
below. 

The philosophy of this approach is that, like in eqtns 11-13, a simple 
Maclaurin series of either 1 or n terms about a Maxwellian is hopeless because 
such a procedure neglects the global distribution of phase density along the tube 
of force. The tails of f(w) and their skew cannot be accurately predicted by a 
one term Maclaurin perturbation expansion of f about a local Maxwellian. It is 
the boundary conditions at the origin and at infinity in the mathematical example 
and at the sun and infinity in the physical example, as well as the strong speed 
dependence of the Coulomb cross section that force a more useful approximation 
scheme. Both from considerations of the conservation equations and from integral 
formulations for the local form of f(w) (SO I, II and Olbert (1983)) global 

correction must be incorporated in the lowest order ansatz for f^°^(w,x) before 
attempting to determine the skew and, therefore, the net heat flux (8b7. In 
mathematical language this represents an attempt to renormalize the perturbation 
expansion of CBES. 

In the face of significant non-local corrections to local gaussian behavior 

we have adopted a lowest order functional ansatz for f°(w, x) that while still 
only speed dependent contains an additional shape parameterT k , beyond the shape 
factor in a gaussian identified with the temperature. This is a significant 
departure from CBES postulate (3). This ansatz is called the kappa distribution: 


( o ) o ~ A n 

f (w.x) = f (k(x),w/(x), n(x), w 2 ) = — — 

(/IH. 


3 ( 1 + w 2 /( k w c 2 )) -(k+1) (33) 


originated by Olbert et al . (1967) and Olbert (1969). This lowest order 
distribution contains extrathermal tails which can support the global Coulomb 
physics, but asymptotes to a gaussian/Maxwellian distribution as c + », In every 
respect the kappa function represents a mathematical "extension" of the Maxwellian 
function. The physical meaning of n is the usual number density, w is the speed 
of the most probably occurring member of the speed distribution, w ?s the proper 
frame speed, and <(x) controls the local strength of the extrathermal tails: small 
k implies strong tails, k going toward infinity implies increasingly weaker tails. 
A k (kOc)) is a normalization constant which ranges between 0.75 to 1. Paralleling 
in concept the Pade approx imant (14) for g by choosing a, b, c and d 
appropriately, <(x) must be determined. This is done self-consistently within the 
perturbation scheme described below. 
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Clearly, as in any perturbation theory, the better the ansatz, the better the 
approximation becomes. In our choice of ansatz we have been influenced by economy 
of free variables, the energy dependence of the observed f( w) , the global aspects 
of the Coulomb physics (as outlined in SO I, II and previous paper) and that a 
member of the functional ansatz should be a gaussian. 

We have used a speed dependent Krook collision term as in SO I and II and 
have successfully imposed conservation of electrons through a mandated k(s). The 
imposition of this eollisional invariant constraint is well known; the equation of 
constraint relates k to other state variables (cf. (34) below). In the presence 
of a speed dependent Krook collision operator this equation of constraint relatfes 
the variation of the non-thermal tail index, d<, to the variation of the 
normalization of phase density d tn (n /w J ) along the tube of force. This 
equation has the form 

d<(s) = - fc(tc(s) ,< (s )) d in (n (s)/w ^(s)). (34a) 

o o e c 

For reasonable interplanetary profiles < is smaller near the sun than at 1 AU. 

This implies that the isotropic part of f(w) near the sun has more of a 
non-thermal tail than at 1 AU. 


Equation (34a) may be recast to show its non-local character 


where v(s) 



v(s)/v(s ) = exp 
o 


:(s) 


-/ 


K (S ) 

o o 


is proportional to n e / w c ' 


die' -i 


(34b) 


K . k q ) ] 

Notice that equation (34b) 


(implicitly) prescribes k(s) as a function of the contrast of collision frequency 
t/hich equals the contrast in the phase space density (at zero proper frame energy) 
between s and some reference point, s , where k= < is known. Mathematically we 
have succeeded in instilling sufficient information in (renormalizing) the lowest 

order f^°^ so that the global, velocity dependent, Coulomb moderated mobility of 
the electrons does not violate charge neutrality. This lowest order failure of 
the CBES, Maclaurin series approach to obtaining extrathermal tails on f(w) is now 
removed . 


It is crucial to appreciate that it is the skewness of the increasingly 
non-local trans and extra-thermals which determine the heat flux (8) for 
astrophysical plasmas. The "resummation" of non-local effects embedded in our 
ansatz (33) as constrained by (34b) heavily impacts the dominant speed, w* , where 
subsequently obtained skews of f make the principal contribution to Q„ (8). This 
effect is exacerbated by the high power. of w involved in determining Q and the 
fact that the perturbation 6f about f itself has significant speed 
dependence . e 

We then obtain a gradient expansion about the local, isotropic f ^ to 
obtain the skew correction term that represents the flow of heat. The analogue of 
(9a) has the global-local (GL) structure (if there exists no current in the proper 
frame which determines the proper frame parallel electric field): 
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Q„ GL (s> 


d d 

-CTS(<( s) , K (s )) - in T - CSSU(s), k (s )) — in (n/T 3/2 ) (35) 

ds ° ds 


CTS > CSS > 0; kT = - mw 2 , 
j* c . c 


where CTS and CSS are global functions of system parameters, which sense, through 
equation 34b, the collision frequency contrast of the physical system connected to 
the observer, while the gradients sense local changes. in the state variables. The 
functional dependence of the global-local form of Q„ (35) should be contrasted 
with the completely local form of CBES in equation 9« 

7 . The CBES Limit of Global-Local 


We recall that a condition via equations (25, 26) for internal consistency of 
the CBES regime in the presence of flows (K * 0) is that 


X — fcn(n/T 3/2 ) = 0(K 2 ) . 
ds 


This implies that the second term of (35) is quadratic while the first is linear 
in the expansion parameter K. [As an aside CTS ^ CSS]. The imposition of near 
adiabaticity implies large collision frequencies, v = 0(1/K), which implies 
progress toward Maxwell-Bolt zmann distribution. Therefore, from the equation of 
constraint (34) since £(k,k o ) is bounded and 0(1) 

< - k(s q ) = 0(K 2 ) , 

which implies that the CBES regime is characterized by a constant 

2 

K S K + 0(K ) 00 , 

O 

which by construction of our ansatz is functionally possible. 

Since lim T z T 
<-> c e 

we may now form the CBES limit of Q m GL (eqtn 11) as 
GL d 

““ Q m = " lim CCTS(k ,k ) - In T + 0(K 2 ) ] = Q„ CBES + 0 ( K 2 ) , 

+ CBES <.-<**> ds 

(K+0) k +» 

o 

where the quadratic terms throughout the CBES approach have already been 
neglected. We reemphasize here the extremely narrow classes of astrophysical 
problems (delimited by equation (29, 30)) for which CBES results are quantitative. 
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8 . The Global-Local Evaluation of the W Function for High Speed Wind Conditions 


As a final result we illustrate in Figure 1 the dimensionless heat function 
W(p ) defined by (16) and central to the previous paper by Olbert (1983) « bub here 
evaluated with Q„ from equation (35). We have used the same coronal hole 
density and magnetic profiles as used in Olbert (1983) for the high speed solar 
wind (U = 700 km/s) and a composite electron temperature profile that has a 
coronal^hole value at 1.03 r (8 x lO^K), has the correct dinT/di,nn 
characteristics near 1 AU (Sittler and Scudder (1980)), and an electron n 
temperature characteristic of the high speed solar wind at 1 AU of 8 x 10 °K. We 
emphasize that these results are preliminary as discussed below, since all aspects 
of the calculation are not fully self-consistent, such as picking a density and 
temperature profile that are actually compatible with the postulate of steady 
state. In that respect the results of Figure 1 represent a feasibility study for 
investigating qualitative regimes of W( p ) that are possible without doing the 
quadruple n-umerical integrals as have been painstakingly done in the previous 
paper. This matter of complete self-consistency is an issue in both these papers 
and is currently under investigation. 

We note the properties of W that the previous paper (Olbert (1983)) reports 
as essential for a conduction supported high speed wind: 1) the initial plateau at 
large densities; 2) the plateau magnitude (this implies sufficient random energy 
available to be converted to directed energy for 700 km/s wind) and 3) the sharp 
decrease near the critical point. The second plateau at lower densities has 
little effect on the dynamics of the solar wind expansion. Figure 2 illustrates 

GL. 

the relative contributions of the two terms of equation (35) to the final Q„ i 
note that different density regimes reflect different relative importance of the 
two terms and neither term is negligible and that the second term of (35) 

(neglected by CBES) determines a larger share of the heat than the conventional 
grad T term at the base of the coronal expansion. 

Figure 1 represents a realization of W(p ) using the hybrid-global-local 
approach instead of the quadruple numerical integrals required in the integral 
formulation of the previous paper. In the preceding paper (Olbert (1983)) a H 
gradations of global Coulomb effects are accounted for in the presence of a 
leading term Ohm's law electric field. In this paper the electric field has been 
set by the zero current condition, particles are conserved and there is a hybrid 
global-local approach to the tail/skew problem. Both approaches attempt, with 
different levels of precision, to assess the importance of the strong velocity 
dependence of the Coulomb cross section which invalidates the CBES approach. 

The W( p ) function is very sensitive to the profiles of n and T g . For a given 
density and magnetic tube geometry the electron temperature profile is not free to 
be prescribed, but is determined by the energy equation (15). In future this 
internal self-consistency will be addressed. For the present these results 
represent significant progress toward affirming the existence of stellar winds 
with large asymptotic velocities supported by Coulomb mediated heat conduction, 
with the solar wind as a prime example. 
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LOG (10) CW = Q/ (*BU*»2) 3 


GLOBAL-LOCAL CHAPMAN-ENSKOG 


LOG (10) DEN 


Profile of the dimensionless heat function W = Q„/(aBU^) for coronal 
hole magnetic field and density profile used by Olbert” ( 1 983) , but 
evaluating Q f( via the global local approach outlined in text. 
Abscissa is electron density. 




REGIMES OF TEMP VS ENTROPY DOMINANCE OF HEAT FLOW 


Fig. 2 
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The fraction of the suggested heat flux determined by the first 
term of the full heat expression of equation 35 for the heat flow 
profile of Figure 1. Note that there exists varying domains of 
relative importance with the second term of (35). not included in 
the CBES expression of equation (9) • playing the dominant role near 
base of coronal expansion. 
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